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Linked cluster expansions provide a useful tool both for analytical and numerical investigations of lattice field 
theories. The expansion parameter is the interaction strength fields at neighboured lattice sites are coupled. 
They result into convergent series for free energies, correlation functions and susceptibilities. The expansions 
have been generalized to field theories at finite temperature and to a finite volume. Detailed information on 
critical behaviour can be extracted from the high order behaviour of the susceptibility series. We outline some of 
the steps by which the 20th order is achieved. 



1. SERIES EXPANSIONS 

Realistic physical models evolve enormous 
complexity and in many instances force us to 
use approximations or expansions about solvable 
cases. Convergent expansions allow for precise 
quantitative statements within their range of con- 
vergence. Under general conditions they are suit- 
able for an investigation of phase transitions and 
of critical phenomena. The price to be paid is 
that normally high orders have to be computed 
in order to reveal collective behaviour, and special 
techniques have to be developed and applied. 

The linked cluster expansion (LCE) amounts 
to a convergent hopping parameter expansion ap- 
plied to the free energy or connected correlation 
functions of lattice models [Q. To be specific, let 
us consider the O(N) symmetric scalar model de- 
fined on the hypercubic lattice A by the action 
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2k , x, y nearest neighbour, 
, otherwise, 
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and lattice site action V($> x ) appropriately boun- 
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ded from below. With the free energy W defined 
by 



we obtain the LCE as the Taylor expansion 



W{v) = exp I E v xy 
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and the obvious generalization to connected cor- 
relations and susceptibilities 
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X 2 ,X 3 , ...,X2n 



The rapidly increasing combinatorical complex- 
ity can be managed by the approrpiate graphi- 
cal device. The application of the LCE in field 
theory has been poineered by Liischer and Weisz 
||. Liischer and Weisz studied a lattice <& 4 - 
theory close to its continuum limit in four dimen- 
sions. They succeeded in a 14th order computa- 
tion of susceptibility series for arbitrary quartic 
couplings. 

Normally, these expansions are applied to mod- 
els on a lattice that is unbounded in all direction, 
and to the study of 2nd order transitions. There 
are recent new fields of application and general- 
ization of LCE. A major topic are field theories at 
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finite temperature |3| . Here one dimension of the 
lattice is compactified to a torus of length given 
by the inverse temperature Lq in lattice units. 
If L is an even number, equivalence classes of 
graphs can be introduced that do not depend on 
Lq and hence are the same as for Lo = oo. There 
are then mainly two adjustments necessary in or- 
der to generalize LCE to this case. First, imbed- 
ding numbers change. The second point is more 
laborious. Comparing field theory models at zero 
and non-zero temperature, the shift of transition 
regions typically is rather small in bare param- 
eters, whereas critical behaviour such as critical 
exponents in many instances undergo a consid- 
erable change. This different critical behaviour 
has to be resolved by the different high order be- 
haviour of the susceptibility series. It amounts to 
an even higher order computation than on lattices 
unbounded in all directions. Graphs contributing 
to the expansion must be able to wind around the 
temperature torus sufficiently often in order to 
pick up sufficient finite temperature dependence. 

We briefly mention two other recent general- 
izations of LCE. The first one is to apply LCE 
in a finite volume and to utilize finite size scaling 
analysis in order to study both 1st and 2nd or- 
der transitions, as well as tricritical behaviour Q]. 
This is the topic of the talk given by H. Meyer- 
Ortmanns, these proceedings. The second gener- 
alization applies to models that couple fields at 
lattice sites which are not nearest neighbours ||. 

2. ALGEBRA OF GRAPHS 

In the range of orders considered, the approxi- 
mated computational costs of increasing the LCE 
series by two orders are a factor of about 30-40. 
Fast algorithms have to be developed to achieve 
this computation. Some of the problems that one 
is faced with are briefly described in the following. 

A first step is to restrict the classes of graphs 
to be considered and write correlation functions 
in terms of quantities that allow for expansions 
into those restricted graph classes, such as 1PI 
ones. Furthermore, 1PI graphs are composed of 
1VI graphs S and renormalized moments Q (one 
external vertex only). 

Incidence matrices provide a convenient way to 



represent graphs algebraically. For a graph T with 
V vertices, the incidence matrix Ir is the symmet- 
ric V x ^-matrix obtained as follows. 

• Enumerate the vertices in any order. 

• Ir(h i) = number of external lines attached 
to the vertex Uj. 

• For i ^ j, Ir(i,j) — number of common 
lines of the vertices v i and Vj 

As an example, the graph 
2 




is represented by the incidence matrix 
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Ir = 



symm. 



V 



This way of representing graphs is not unique be- 
cause the vertices can be enumerated arbitrar- 
ily. Comparison of two graphs needs huge fac- 
torial numbers of simultaneous permutations of 
rows and columns. The actual construction of 
graphs is done by iteration on the number of lines. 
Equivalent graphs have to be kept only once, and 
the number of comparisions is a measure of the 
efficiency of the algorithms. 

A " canonical" representation of a graph T that 
is both unique and efficient is the following. In- 
troduce a (partial) order relation on the set of 
vertices and enumerate the vertices in accordance 
with it. The incidence matrix then is defined by 



7^"(i,j) = max/ r (7r(i),7r(i)), 

7rEil' 



(6) 



where II' is the subset of permutations that ex- 
change vertices only if they stay unordered rel- 
ative to each other under the preordering. The 
maximum is with respect to any total order rela- 
tion on symmetric matrices. 

The preorder of vertices should be as complete 
as possible. If all "local" properties of vertices 
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are taken into account in an optimal way we ob- 
tain the order relation introduced by Liischer and 
Weisz. They succeeded in a 14th order computa- 
tion for 2- and 4-point functions. Even higher or- 
ders need to take into account more global prop- 
erties. A good example is given by 

3 3 3 3 3 




with many 2-vertices staying unordered. Taking 
into account properties like the distance to ver- 
tices that have more than 2 lines attached, vertex 
ordering can be enhanced to give 

4 6 7 5 s 



2, 


' 4 6 7 5 3N 






V ^l6 7 5y 





There are (3!) 5 permutations left to build I™" 
instead of (15!). Furthermore, order relations can 
be iteratively improved by taking into account 
order numbers of neighboured vertices, next to 
nearest neighbours and so one. Normally, one or 
two iterations are sufficient. 

The table shows the number of graphs for some 
classes generated as yet. The index k = 2,4,6 
denotes the number of external lines attached, L 
is the order (number of internal lines). 

The contribution of a graph, its weight, con- 
sists of two factors. The first factor is a rational 
number, itself a product of internal and (inverse) 
topological symmetry numbers and of the lattice 
imbedding number. It is only the latter that de- 
pends on the lattice size. The second factor comes 
from the product over the vertices. Every ver- 
tex contributes a finite dimensional integral that 
depends on the number of lines entering the ver- 
tex. The complete weight is a real number. The 
accumulation of roundoff errors is avoided by in- 
troducing so-called vertex structures. A vertex 
structure of a graph T is the multitupel of non- 



negative integers, T — > (/ii, jjb%, . . .), where jii de- 
notes the number of vertices of T that have pre- 
cisely i lines attached. Before multiplying the real 
coupling constants of the vertices, graphs with 
identical vertex structure are summed according 
to their rational weights. 

An application to finite temperature field the- 
ory can be found in nj. 



L 


Q,2(L) 


S 2 (L) 


SAL) 


S 6 {L) 


n 
u 


1 
1 


i 
i 


i 
i 


i 
i 


i 
i 


n 

u 


n 
u 


u 


u 


2 


1 


() 


1 


1 


3 





1 


1 


2 


4 


4 





4 


6 


5 





2 


4 


11 


6 


15 


3 


20 


4G 


7 





8 


27 


91 


8 


79 


9 


117 


349 


9 





40 


214 


837 


10 


439 


68 


815 


3140 


11 





247 


1830 


8401 


12 


2877 


470 


6721 


31187 


13 





1779 


17028 


90599 


14 


20507 


3937 


61653 


336582 


15 





14801 


170923 


1042392 


16 


161459 


35509 


621191 


3895341 


17 





135988 


1834324 




18 


1376794 


350614 


5548427 




19 





1361878 


20967387 




20 


12693105 


3705467 







REFERENCES 

1. M. Wortis, "Linked cluster expansion", in 
Phase transition and critical phenomena, 
vol.3, eds. C. Domb and M.S. Green, Aca- 
demic Press, London, 1974. 

2. M. Liischer, P. Weisz, Nucl. Phys. B300- 
[FS22] (1988) 325. 

3. T. Reisz, Nucl. Phys. B450 (1995) 569; 

4. H. Meyer-Ortmanns, T. Reisz, "Critical Phe- 
nomena with Linked Cluster Expansions in a 
Finite Volume" , [hcp-lat 960400(1 

5. A. Pordt, T. Reisz, "Linked Cluster Ex- 
pansions beyond nearest neighbour interac- 



tions: Convergence and Graph Classes" , [hcp- 
lat 96040211 



4 



6. T. Reisz, Phys. Lett. B360 (1995) 77. 



